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Abstract 

A layered material structure in a monochromatic neutrino 
beam produces interference effects that could be used for 
the measurement of features of the neutrino mass matrix. 
The phenomenon would be most useful at high energies. 
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The MSW effect [1,2] describes how electron neutrinos in matter propagate 
differently from other neutrinos, and from electron neutrinos in vacuum. This effect is an 
element in the interpretation of recent experiments [3-6] that have explored the neutrino 
mass spectrum. The phenomenon also describes the effects of the presence of boundaries 
between different media on neutrino propagation and oscillation. As recent work [7-9] 
has shown, the boundaries introduce the possibility of interference between different 
amplitudes for neutrino propagation. Indeed, the presence of boundaries within the earth 
has implications [7,8] for the interpretation of the data of reference 3. 

Here we point out that an arrangement of layers of materials, containing many 
boundaries, can provide another angle on the interference between the propagation 
modes. In addition to describing the basic mechanism for a two-family neutrino structure, 
we briefly address the question of where this effect could be most profitably employed. 



I. Review of MSW and Boundary Interference effects 

Neutrino oscillation occurs because the electroweak eigenstates are not the mass 
eigenstates. Let us consider two families of neutrinos, with electroweak labels e and 
(generically Greek letters) and mass eigenstate labels 1 and 2, with 1 labeling the lightest 
neutrino. The matrix U (the leptonic CKM matrix) connects these states according to 

K>=Ic/„,K.>, (1.1) 



(=1,2 



where U takes the generic form 

COS0 sin0 
-sin0 COS0 



(1.2) 



with Q such that cos(20) is positive. 

Starting with a pure beam of, say, v^, with definite momentum (assumed 
throughout), the time evolution is governed by the mass eigenstates and gives after the 
neutrino beam has traveled a distance x = r a beam containing a mixture of each type of 
neutrino in the usual fashion. In particular the probability for the conversion )S — > a is 

/'MO = sin^20 sinV, (1-3) 
where the angle <p is determined by the energy difference, 

1 Am^ 

(p = -{E,-E,)x = X, (1.4) 

^ 2^ ' AE 

2 2 2 

with Am = m2 -m.\ . The approximation refers to the limit in which the mass difference 
Am is much less than the momentum of the beam, a limit that will be of interest to us. 
The oscillation length €(£) is the distance that corresponds to a change in ^by n\ that is 
€(£■) = AitEllSm. Finally note that the probability for nonconversion, i.e. that ;S — > ;S, is 

Ppp = cos^^ + cos^20 sin^9 = 1 - P^a- ( 1 -5) 

Propagation in matter. In the presence of matter, each neutrino specie may have a 
different potential. (In normal matter it is only the electron neutrinos for which there is a 
potential associated with rescattering from electrons in the material.) In particular, 
suppose that the potential difference of the two neutrino species is Vk, where k labels the 
material. Then the effect of the matter on the propagation parameters is described by 
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where 



Am^ Ami = Am^ ^(cos 29 -^^f + sin^ 29 



Tk k 



9^9,^ sin' 29, = 



(1.6) 



sin' 26 



(cos 261 -^,)+ sin' 29 



Am' 



(1.7) 



The value of E for which = cos20 is referred to as the MSW resonance; for that value 
the angle 20, goes through 71/2. for general values of 0, total conversion in a single 
thickness of any medium, including vacuum, is not possible, but it can occur for an 
appropriate thickness of medium at MSW resonance. 

The expression for 0, above is equivalent to the form given in Ref. [2], namely 

tan 20 „^ 

""^■ ^ 4.E 

1 ^ sec 20 

^oAm' 

with the replacement €o = 27t/Vi. 

Let us use the label k for the material (carrying information not only on the 
composition through the neutrino potentials but on the layer thickness as well). Then the 
generic amplitude ApJ'^^ that a p-neutrino enters and a (7-neutrino leaves layer {k} is 
given for the two neutrino types a and P by 



: cos (p, + i cos 29, sin (p, , 



A^p^ = cos (p, - i cos 29, sin (p, , 

4? = <^ = -''sin26l,sin9,. 

Since we shall mainly be concerned with total conversion it is useful for later 
comparison to give here the length Xk for maximum conversion in a single layer of 
material k, immediately found from the last of these equations. The maximum possible 
conversion is realized for a length such that the factor sin^^ = 1, i.e.. 



(1.9) 



X,= 



4E 7t 



Am, 2 



Am K 



Am, V, 



4 



7t 



cot 20 



7t 



n 



far below MSW 



at MSW 



far above MSW 



(1.10) 



2 2 

The maximum conversion probability is sin 20^:, which is sin 20 far below MSW, unity at 
MSW, and sin^20/(^^^, asymptotically small, far above MSW. 
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The smallest value of energy and hence of leads to the smallest length X^. At 

the same time, it is easiest to detect energetic neutrinos. It is therefore helpful to have 
some idea of the energies that are involved for neutrinos. We note that when we use the 
number = 6 x 10"^ cm appropriate for earth [2], then 

= (E/(AmcY) X 2.5 x 10"^^ eV 

For Amc^ = O(10"^) eV, is on the order of 1 for £ = 0(4 MeV). One can think of this 
energy as roughly the dividing line for whether one is below or above MSW, although of 
course the precise MSW energy depends not only on but on the correct value of Am as 
well, and the latter number is not yet fully understood [10]. 

Propagation through layers. When there are layers of matter with differing 
densities, then interference is possible. The history of the effects of passage through 
repeated layers on mixing is in fact extensive, beginning perhaps with a discussion of 
neutron- antineutron oscillations in a nonuniform magnetic field [11], then continuing 
with work on neutrinos by Ermilova et al. [12] and later by Akhmedov [13] in a long 
series of papers. These latter references describe what the phenomenon of parametric 
amplification; in fact the work of refs. [7-9] can be seen as a special case of this treatment 
[14]. However, the approaches taken in refs. [7-9] and in this paper are oriented in a 
fashion sufficiently different to make them worth independent consideration. 

Let us consider two layers oriented perpendicular to a neutrino beam, the first 
labeled { 1 } and the second {2}. Then, as pointed out in [7], the amplitude Aap for passage 
through two successive layers { 1 } and {2} (what we refer to as a bilayer) contains two 
terms, and these terms can interfere: 

= sin(202~20i)sin^iSin^2 (l-H) 

-zjsin 20, sin^j cos^j +sin 202 cos^, sin ^2} 

One can see immediately that the structure of this amplitude is not that of the single layer; 
for example, the third of Eqs. (1.9) is purely imaginary. In particular it is possible for this 
amplitude to have magnitude one — total conversion — over a wide range of the parameter 
space. 

The work of reference [7] approaches total conversion for a two-channel problem, 
with a beam initially of type a passing through a double layer, through the probability 
condition IA„^P = 1. In a rather involved calculation it is shown using this condition that 
total conversion occurs for layer thicknesses such that 

2 2 -cos 20, 



cos2ft cos(202-2a) 

' (1.12) 

2 7 — COS Zfi 

^2 = tan^ ©2 = 7 — 7 

COS 202 COS (26>2- 2611) 

in regions of 0i and 02 where the right hand sides are positive. (Note our definition = 
tan(pi.) 
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In the two-channel problem Eq. (1.12) is in fact more simply approached through 
the amplitude condition A = 0. This amplitude is given by 

= COS 9i COS (P2 - COS {262 - ) sin sin (1-13) 
-/ {cos 20J sin 9j cos (P2 + cos 202 cos (p^ sin }• 

This amplitude is zero when both its real and imaginary parts vanish, representing two 
conditions for the two angles (pi and ^ and hence for the layer thicknesses if all other 
physical parameters are given. One can see immediately from Eq. (1.13) that these 
conditions are simply written as conditions for yi and y2, namely 

Real part = 0: 1 - 3;i3;2cos(202 - 20i) = 0. (1.14a) 

Imaginary part = 0: };icos20i + 3;2COs2^ = 0. ( 1 . 14b) 

From these equations, quadratic in the y,, one immediately arrives at the solutions given 
inEqs. (1.12). 

II. Passage through multiple layers 

The amplitude Aaa for the survival of neutrino type a through multiple layers is 
developed in a straightforward way from the single-layer amplitudes of Eq. (1.9). In the 
two-family problem, this amplitude is an element of a 2 x 2 matrix resulting from the 
multiplication of two primitive (single layer) 2x2 matrices. The generalization to more 
than two layers is straightforward. We give here the cases of three and four layers as 
examples; in each case we give the conditions for = 0. 

For three layers, the conditions that the real and imaginary parts of Aaa = are, 
respectively, 

l-t,yiyjcos{2d,-2dj) = (2.1a) 



cos 20,. - y^y2yT, cos [26^ - 2^2 + 2^3 ) = . (2. lb) 

For four layers the respective conditions are 

4 

1 - ytyj cos (20. - 2d. ) + y,y2y3y4 cos {2d, - 262 + 2d^ - 2^4 ) = (2.2a) 



i.M 
i<j 

4 



J^y,cos20,- Y >',y,y,cos(20,-20,+20,) = O. (2.2b) 



i<j<k 



These two examples are sufficient to understand the more general cases. The only 
important feature to note here is that for more than two layers the two conditions that A«« 
vanish are insufficient to determine uniquely the yt and hence the layer thicknesses. 
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Repeated layers. A solvable case is that of alternating layers with every other 
layer identical to its partners. In other words, we have exactly repeating layer pairs or 
repeating layer pairs plus a last layer identical to the first. If we label A'^ as the total 
number of layers, then these possibilities correspond to even and odd, respectively. 
The number of bilayers is n = [Nll'\, where the square bracket indicates the largest integer 
in Nil. We found no especially interesting solutions for the odd A'^ case and will make 
only passing comments on it. 

For even A^, we have in mind ultimately a situation in which the first member of a 
bilayer is vacuum and the second is a given thickness of a dense material, but we treat the 
more general situation of a separate potential difference for each layer. In this case Ami = 

2 2 2 2 

Am3 = Ams = . . . ; Ami = Am4 =. . . ; = ^ =. . . ; and so forth, so that we have only the 
subscripts 1 and 2. The two conditions for real and the imaginary part will now determine 
yi and y2. 

We give a series of explicit results for the conditions for total conversion for 
multiple bilayers in the appendix. We remark here that the 2 x 2 matrix A that gives the 
amplitude for the passage through n bilayers can be written as a factor (cos^i cos^)" 
times a remaining matrix A'. Since the conditions refer to the vanishing only of the aa- 
component of A, we derive the (necessary and sufficient) conditions from A'^a = 0. These 
are the conditions given in the appendix. 

The calculations presented in the appendix reveal two important feature that we 
shall assume to be general: First, the imaginary part of the amplitude A'aa contains a 
single factor of the combination 

F = y\co^29i + j2cos202 (2.3) 
This will turn out to be quite useful, as we shall see below. Second, aside from this single 
factor, the modified mixing angles 0, appear in A'aa only in the combination 5 defined by 

5=201-202. (2.4) 

Conditions for total conversion. In the case of the single bilayer, the imaginary 
part in particular vanishes only if the factor F defined by Eq. (2.3) vanishes. For more 
than one bilayer, either F or its coefficient could vanish. Let us consider the latter 
possibility for some low order examples. 

For n = 2, the imaginary part of the amplitude is given by Eq. (A. 2b), and we want 
to consider the possibility that the second term vanishes, i.e. that 1 - yiy2C0^5 = 0. This 
gives y2 - (yiCOsd)~\ and when this is substituted into the real part, Eq. (A.2a), we find 
the condition 

COS O COS o 

But each term on the right side of this expression is positive, so we do not have a 
solution. 

For n = 3, we consider the possibility that imaginary part vanishes because the 
expression in curly brackets of Eq. (A.3b) vanishes. But three times the curly bracket in 
Eq. (A.3b) + the left side of Eq.(A.3a) (the real part for n = 3) is 

-^Ui +^)—2 — rr- 
cos o 
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and, as for its analog in n = 2, this quantity cannot vanish. 

Although once again we do not have a general proof, it is reasonable that the only 
way for the imaginary part of the amplitude A'aa to be zero is with the condition that F 
vanishes. Using this condition we can make an arbitrary n generalization for the form of 
the amplitude at the total conversion point. To do so we write the amplitude for passage 
through a single bilayer in canonical form, namely 



B = 



e'""" cost; e"""" sinr; 
e ^ smTj e cosrj 



(2.5) 



The parameters of this "unit cell" amplitude are determined by comparison to the explicit 
result, 

where the single layer amplitudes A'^' are given by Eq. (1.9). Using Eq. (1.9), we find 
that B has the more restrictive form 



B = 



(2.6) 



e COST] e ^ smrj 
-e'""^ sin 77 e cosr] 

with the three bilayer parameters of this expression, r\, cOaa, and cOap, given in terms of the 
single layer parameters by 



cos 77 



1 



(1 - COS d) = cos cos ^'2 (1 ~ ^2 cos 5) (2.7) 



tan«„3 =- 



cos 20^ + cos 202 
1~ JiJi cos 5 
sin 20J + ^2 sin 2^2 



(2.8) 



(2.9) 



y^y2 sin 5 

This expression simplifies further if we apply the condition that for n bilayers the 
factor of Eq. (2.3) is zero at the total conversion point: 

O)aa = for total conversion (2.10) 
(The denominator of Eq. (2.8) is not independently zero except for the single bilayer 
case.) 

With this condition, the n bilayer amplitude with total conversion of the a beam 
becomes 



B" = 



cos nrj e ^ smnr; 
-e ^ smrO] cos nT] 



In turn, we see immediately that for total conversion 

cos(nri) = 0. 

In turn Eq. (2.12) gives 

(2m + l);r 



7]. 



2n 



m = 0,l,...n — l 



(2.11) 

(2.12) 
(2.13) 



The pair of conditions that F and cos(nr7) each vanish provide us with two 
equations for the angles (pi and (jh, i.e., for the layer widths xi and X2. We shall describe 
the solution to these equations in the next section. 
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III. Total conversion in a repeated multilayer system 



We apply the simultaneous conditions F-0 and cosnr] -Oto detennine yi = tm(pi here. 
The F = condition determines yi in terms of yi- The quantity yz is determined in terms 
of the angle rj by the inversion of Eq. (2.7), which is a quadratic equation for y2 in terms 
of COST]. Since cosnj] is an nth order polynomial in cost], there are 2n solutions for 3;2- 
Because A'^ = 2n, this matches the number of solutions coming from the Mh order 
polynomial for y2 coming from the original real part equations, as described below Eqs. 
(2.10). Thus we find all the solutions in this way. Because we would like to minimize the 
thickness of the material layers, we shall be interested in small solutions, and we shall 
see that this corresponds to small values of rj. 
If we define 

Zi = yiCos2di, (3.1) 
then the condition that F vanish reads Zi + Zi = 0. In turn this means that 

Zi^ = Z2^^z^. (3.2) 

Equation (2.7) now becomes 

cosT] = ^ ^ (cos20icos202 + z^cosg) (3.3) 

^(z' + cos' 2d, ) (z' + cos' 202 ) 

One can see quickly that for small J], for which cosr; 1, Eq. (3.3) becomes 

2 2 

homogeneous in z , and so has solutions at z =0. For more detail, we consider separately 
different regimes of the MSW parameter <^,. In doing so it is simplest to treat our the layer 
labeled 1 as a layer of vacuum (^i = 0). It is straightforward to generalize to a bilayer 
consisting of two different materials each with nonzero values of ^i. 

^2 small (below MSW resonance). We treat ^2 as a small perturbation, with Ami^ 
= Am^ and 62 = 6. We have 

cos(5= 1 -(1/2)^2 W20, 
and hence Eq. (3.3) becomes to leading order in ^2 

(1+3^2 cos' 26*) tan' 20 ^ 



/ ^l-cosr; = 



2(1+ ^2)' 



^2 

With the definition a = -^j , this equation reads 

yl (1 + yl cos' 20 )a tan' 20 = (l + y] )' . (3.4) 

The right side of this quadratic equation for ^2 is 0(1) or larger, so there are no solutions 
unless a = 0(1), i.e.,/= 0(^2^). But by comparison with Eq. (2.13), we see that for small 
values of m/n, 

1 (2m + lf7t^ 

f = -- i . (3.5) 

2 4n 

By choosing n large enough, or more particularly m/n small enough, we can imagine 
choosing a = 0(1). The equation for y2 will then have small positive solutions. 
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The formal solution of Eq. (3.4) is 



^ ± sin^ 26 ^a(a-4 cos^ W)-2 cos^ W+a sin^ W 

" 2 cos' 20(1 -a sin' 261) ' ^^'^^ 

We see immediately that there is no (real) solution unless a > 4cos 29. The suitable (both 
small and positive) solutions correspond to the minus branch y2^-, and it is this solution 
that we look at henceforth. The solutions are smooth as we pass through the point where 
the denominator factor 1 - sin 20 = and are simple for a greater than or equal to 0(1) 
for the 6 = 0.7 case that we use below for illustration. Indeed, in this range we can use the 
very accurate approximation 



2 cot 26 2^^\ 

y\_= + (2-cos' 26) 

a ^ ' 



c ot' 26 
a 



+ . 



2 

all the way to small values of a. That is because the expansion is in cot 261a, something 
that follows from the fact that it is Otan 26 that appears in the original equation. Thus we 
can write our solution in the form 

_ cot 20 _^ 'Kco\.26 
4a 2n^2 

where in the last step we have chosen a large n solution with m = 0. Under the assumption 
that cot20 is small enough, we can replace by (jh, and solve for X2: 

X. = cot 20 . 

The total amount of material nx2 is less than the corresponding amount of a single layer of 
material 2 [Eq. (1.10)] only in the circumstance that cot20is very small. 

There is a second problem in this region of ^2- Once we have found y2, then yi, 
and hence the thickness of layer 1, is determined through the condition 

= JiCOs20i + J2COS202 

But in this region 61 = 62= 6. We would then require yi = -y2, or, assuming that y2 is 
sufficiently small that 92 = tan<3!>2 = y2, 

(pi = tan"^(-92) 

Since the distances X2 must be positive, the only way we can satisfy this condition is to 
take (pi = 27t- 92, or in other words, 

Xi = -—^-X2. (3.7) 
Am 

The first term on the right side of this expression, which is the oscillation length in 
vacuum, is not necessarily small. The multiple bilayer arrangement offers no advantages 
below MSW. 

^2 at MCW resonance. We have 61 = 6 and (pi = (Am^/4E)xi. For medium 2, <§2 = 
cos20, and sin202 = 1 or 202 = ^2 and cos202 = 0. (In fact we shall assume that we are a 
little above MSW resonance, so that cos202 is small and negative. This helps to clarify 
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limits.) The angle ffl, = — ^;c, = x. . We also nave o - ml - Id and cosd = 

sin20. In this limit the relation between y2 and 1] of Eq. (3.3) becomes 

cosr7= J— (3.8) 

The relevant solution to Eq. (3.8) is 

yi = tant]. (3.9) 

or, since y2 = tan^, 

n(2m + l) 

92 =^ = -4 -■ (3.10) 

2n 

Since cos202 is small and negative, yi is satisfactorily positive and is also small: 

j'l = j2lcos2e2l/cos2a (3.11) 
When we calculate the layer thicknesses we see why this case is of no special 
interest. Recall [Eq. (1.10)] that at MSW resonance a single layer of thickness X2 - (jtlVi) 
cot20of material 2 gives total conversion. In comparison, Eq. (3.10) shows us that the 
minimum value of (fh for an n-bilayer system occurs for m = 0, in which case we have (jh 
= Kl{ln), or 

It AE \ Tt i-,x 
^ = cot20. (3.13) 

2n Amj n V2 

Thus the total amount nx2 of material 2 is exactly the amount needed for the single layer. 
Moreover in the MSW limit, cos202 is zero, so that from Eq. (3. 1 1) y\ and hence the total 
thickness of vacuum nxi vanishes. The entire system limits to a single layer of material 2. 

We have been able to find no quantity associated with n bilayers that scales to any 
experimental advantage in the MSW limit. 

^2 large (above MSW resonance). Again we assume that the medium labeled 1 is 
vacuum, ^1 = 0. For ^2 » cosSand sin0, sin^202 = ?>ir^^26l^2 0, with sin202 positive 
and cos202 approximately -1. The fact that sin202 is small means [Eq. (1.9)] that one can 
at best have very little conversion in a single layer of material 2. Thus the very possibility 
of total conversion makes this limit interesting. 

The relation (3.3) between _y2 and r\ gives solutions independent of <§2 in the large 
^2 limit, namely 

sin cos 20 



J2± 



= +- 



^|cos^ 20-cos^r7| 
The positive y2 solution is then simply 

Isinrjlcos 26 

y2= ,| ' ' , , (3.13) 

J cos 29 - cos 77 

Before we deal with the issue of many bilayers, let us consider the case of a single 
bilayer. We show here that total conversion may not be possible in the single bilayer, 
although it will always be possible for n>2. The original total conversion conditions for 
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the single bilayer are given by Eqs. (1.14). If we take yi from the second of these 
equations and substitute into the first, we find an equation for yi, namely 

, cos 20, cos (20 -20,) 

l+y' - = (3.14) 

cos 26/ 

Well above MSW, 202 = ^ - e, and expanding to leading order in e gives 

yl= ^- • (3.15) 

' l-etan20 

For this equation to have a valid (positive) solution, one requires that e tan20 > 1, and this 
will not always hold; indeed it can hold only in a decreasing domain of as £■ becomes 
larger (e 0). This situation is illustrated in the numerical example of the next section. It 
is not difficult to show that there will always be a total conversion solution in this limit 
for two or more bilayers. 

Let us turn next to the case of many bilayers. Given that we are interested in the 
case of small t], and supposing that cos20 is much larger than sinr;, we can replace the 
denominator in this expression by sin20, and our solution becomes 

^2 = cot20 sin77. (3.16) 
In turn, Eq. (3. 1 1) gives us the (small) value of yi, namely yi = J2ICOS202I/COS20 = 
sinr7/sin20. (In the numerical example treated in the next section, we chose = 0.7, in 
which case cot20= 0.17, while sin20= 0.98.) 

We have Ami^ = Arn' ^2 in this limit. We also choose the minimum value 7tl{ln) 
for t], and expand for small 77. Then we compute from our results for yi and y2 the total 
amounts of material 2 and of vacuum space to be, respectively, 

K Tt 

nxr,= — cot 20 and nx.= — csc20. (3.17) 

These numbers should be compared to the length of material X needed for maximum 
conversion in the large ^2 limit, namely [Eq. (1.10)] X = 7t/V2. We see that if the amount 
of material is the controlling issue one can gain considerably, in that one may have nx2 « 
X. However, nxi » X, so that if the total length of the experiment is the controlling issue 
this limit is not useful. We should also recall that the maximum conversion in a single 
layer of width X is asymptotically, so the very possibility of total conversion is an 
attractive feature of the multiple bilayer arrangement. 



IV. Numerical example 



As indicated by the discussion of the previous section, the most interesting cases to look 
at are those for which ^2 puts one above the MSW resonance. We present two numerical 
illustrations here, each for the arbitrary value of = 0.7, corresponding to a large degree 
of mixing. Our bilayer consists of a layer of vacuum followed by a layer of a material 2 
for which the potential is given by V2 = 6 X 10"^ cm [2]. In the first example we assume 
the energy and masses are such that one is slightly above MSW and in the second 
example one is far above MSW. Our strategy is to first allow the possibility of total 
conversion by fixing the thickness x\ of the first layer in terms of the second layer through 
the condition that F, as defined in Eq. (2.3), vanish, i.e. through Eq. (3.1 1). We then plot 
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the probability for nonconversion as a function of the total width X2 = nx2 of the material 
layers for various numbers n of bilayers, including the single bilayer. 

For the first example, we suppose that we are slightly above MSW resonance, 202 
= kI2 + 0.02. Figure 1 shows the total length of material versus the probability of 
nonconversion (i.e., total conversion is a zero in this plot) for n = 1, 2, 3, and 6. There is 
very little dependence on the number of layers. For the parameters used one can directly 
locate the first large-n zero [Eq. (3.13)], and it matches the numerical value on the plot 
precisely. We have also made a variation on this calculation, in which we have shifted x\ 
from the zero-determining value, and we have observed the zero fill in as the shift 
increases, verifying that the conversion is no longer total. 

As a second example we suppose that we are far above MSW, 202 = K - 0.02. 
Figure 2 again shows the probability of nonconversion for n= 1, 2, 3, and 6 as a function 
of the total amount of material used. In this case the single bilayer does not give total 
conversion. One can see the zero move to the left (less material) as n increases, with the 
overall pattern quite distinctly dependent on n. 

V. Comments 

We have concentrated here on the possibility of total conversion of neutrinos in 
multilayer systems. It would appear that the technique is more interesting at high 

energies. If these effects are ever to play a role in experiments it will be important to 
understand several features that we have not looked at, including in particular the 
implications of a realistic energy spread and, less importantly, the generalization to three 
families. The three family calculation in principle has a richer variety of possible 
outcomes for conversion experiments. 

The neutral ^-system presents another well-known case of oscillation. It differs 
radically from the neutrino system; among other differences materials in the kaon beam 
produce absorption as well as forward scattering. This system may be interesting to think 
about from the point of view taken here. 

Finally we remark that there is another class of effects that exploits the fact that 
the order of layers matters in conversion probabilities. We shall discuss this elsewhere. 
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Figure Captions 

Figure 1. The total length of material, in units of 10 cm, in a multiple bilayer system 
consisting of n alternating slices of vacuum and material versus the probability of 
nonconversion. Total conversion is a zero in this plot. The material has the density of the 
earth, and the relative width of the layers is determined so that the factor F of Eq. (2.3) is 
zero, which guarantees the possibility of complete conversion. The energy of the neutrino 
beam is such that we are slightly above MSW resonance, 202 - Ttll + 0.02. Plots are 
drawn forn = 1, 2, 3, and 6. 

Figure 2. The total length of material, in units of 10^ cm, in a multiple bilayer system 
consisting of n alternating slices of vacuum and material versus the probability of 
nonconversion. Total conversion is a zero in this plot. The material has the density of the 
earth, and the relative width of the layers is determined so that the factor F of Eq. (2.3) is 
zero, which guarantees the possibility of complete conversion for n > 1. The energy of the 
neutrino beam is such that we are far above MSW resonance, 202 = k- 0.02. Plots are 
drawn for n = 1, 2, 3, and 6. 



Appendix 

Here we work through a series of cases of total conversion in n bilayers in order to 
develop insight to the most general case. For n = \, Eqs. (1.14) apply, although it is useful 
to repeat them here. Through n = 4 we find the conditions for total conversion (the "a" 
and "b" equations refer respectively to the real and imaginary parts) 
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n=l: 



1 -yiy2Cos5= 



(A.la) 



n = 2: 



n = 3: 



3;icos20i + 3;2cos202 = 0. (A. lb) 

l-yi^ - y2 - 4_yi_y2COs5 + _yi^_y2^cos(25) = (A.2a) 

l\y\co%lQ\ + 3;2cos202][l - y\yico%5\ = (A.2b) 

1 - 3yi^ - 3y2^ + 3yi^y2^ - 3yiy2Cosd [3 - yi^ - y2^] 

+ 6ji Vcos(2^ - Vcos(3^ = (A.3a) 

[3;icos20i + 3;2cos202] [3 -yi^- y2 + y\yi - 83;i3;2Cos5 

+ 2^1 Vcos(25)} = (A.3b) 



n = 4: 



1 ^2 ^-2, 4, 4,1^22 o2 2/2, 2s, 

1 - Oji - by2 + >'i + J2 + loji y2 - 2yi y2 (yi + y2 ) 
- Syiy2COs5 [2 - 2yi^ - 2y^ + 71^2^] + Vcos(25) [5 - y^ - yh 

- 8 yi Vcos(35) + yi Vcos(45) = (A.4a) 

-2[);icos20i + j2cos202] {-2 + 2);i^ + 2)^2^ - ^y\yi 
+ jij2COs5 [10 - 2yx - 2y2 + ji - 6ji Vcos(2^ + Vcos(35)} = 0. (A.4b) 

We can also write systematically pieces of tenns in X aa for general n. As 
examples, the terms in the real part that are proportional to y^ (by 'f we mean in general 
yi'y-f''', q positive) take the form 

(jij2)"cos(n^, (A.5) 
while the terms in the real part that are proportional to y^" ^ are 



sin(n(5) j^ + y^ sin[(n-l)5] 



yxyi . c + 



sin 5 2 sin 5 



(A.6) 



The terms proportional to y^ ^ in the imaginary part are 



{yiy2T\y\cos2di + j2cos202] (A.7) 
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These terms are the largest powers of y possible in both the real and imaginary parts. The 
real part contains even powers only, with the largest power y^; the imaginary part 
contains odd powers only, with the largest power \ 

Finally, we can also develop systematically low powers of y in the imaginary and 
real parts for arbitrary (large) n. A few examples are 



Real part, constant terms: 
Real part, y^ terms: 



1 



{nll)[lnyiy2COsd + (n - \){yi^ + j2^)] 



Real part, y terms: 



n(n-l) 
4! 



2n {n + 1) ylyl cos 25 +4n{n- 2) y^jj {ji + yl)^^^ 5 
+ {n-2)\{n-3){yf + ylf + 2{n + 3) 



2 2 

yi yi 



(A.8) 
(A.9) 



(A. 10) 



Imaginary part, y^ terms: 

n\y\cos20i + 3;2cos202] (A.l 1) 

Imaginary part, y^ terms: 

"^"^ [ji cos 201 + y^ cos 202 ] [2 (n + 1) 3;i3;2 cos 5 + {n-2)[yl + yl )] (A. 12) 

o 

We have worked out such terms all the way through the y terms. We have not, however, 
found a way to generalize every term. 
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